Two-Hole and Four-Hole Bound States in a t — J Ladder at half-filling 
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The two-hole excitation spectrum of the t — J ladder at half-filling is studied using linked-cluster 
series expansion methods. A rich spectrum of bound states emerges, particularly at small t/J. Their 
dispersion relations and coherence lengths are computed, along with the threshold behaviour as the 
bound states merge into the continuum. A class of 4-hole bound states is also studied, leading to 
the conclusion that phase separation occurs for t/J < 0.5, in agreement with other studies. 



I. INTRODUCTION 

Ladder models of magnetic and electronic systems, 
such as the Heisenberg, t — J and Hubbard ladders, have 
attracted great interest in recent yearsBi, on both the- 
oretical and experimental grounds. The ladders provide 
a "half-way house" between one-dimensional and two- 
dimensional systems, and display new and interesting be- 
haviour in their own right. The spin-| Heisenberg ladder, 
for instance, forms a magnetically disordered spin liquid, 
with a finite spin gapS. Experimentally, two-leg S = h 
ladders are found in some cuprates like SrC^Oai; while 
the ladder material Sri4_ a; Ca a: Cu2404 has recently been 
found to show superconductivity under high pressureS. 

In this paper we study two-hole bound states which oc- 
cur in the t — J ladder at half-filling, using a newly devel- 
oped linked-cluster series expansion method for 2-particle 
statesi. A rich spectrum of multiparticle bound states 
will be revealed. The properties of 2-particle states, such 
as their dispersion relations, coherence lengths and exclu- 
sive structure factors, can provide much important infor- 



mation about the dynamics of the model at hand, provid- 
ing a stringent test of any analytical model of the system. 
These properties are also amenable to experimental test, 
via neutron scattering or other techniques. In the case of 
the Heisenberg ladder, for instance, a modelBB based on 
triplet excitations of a ground state consisting of singlet 
dimers on each rung of the laddcrBSS provides a reason- 
able qualitative description of the numerical results for 
both 1- and 2-particle excitations, not only in the dimer 
limit, but even for the isotropic case. Ideally, one would 
like to achieve a similar convergence between analytical 
theory and numerical 'experiment' for the case of the t—J 
ladder. 



The basic properties of the t — J ladder have been 
explored by exact diagonalizationflli0'0Ei0O, the 
density renormalization group@0@0, quantum Monte 
Carlo simulationil, and series expansion techniquesil0, 
as well as approximate analytic theories^. 
The Hamiltonian of the t — J ladder is 



H — J ^(S i>a • Sj + i, a — -^i,a«i+l,o) + J± ^~)(Sf t i • Si,2 — J«i,l^i,2) 
-t E ^(4a, CT C l+ l,a, CT + H.C.)P ~t ± Y, P(4,l,*Ci,2,a + H.C.)P 



(1) 



where i labels sites along each chain, a (=t or j) and a (=1,2) are spin and leg indices, and P is a projection 
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operator which excludes doubly occupied sites. The con- 
stants J, t are exchange and hopping parameters on each 
chain, while J±, t± are coupling parameters between the 
two chains, i.e. on the rungs of the ladder. 

The scenario is then as followsEl. 

Starting from the half-filled case, and in the dimer limit 
of strong interchain coupling J±,t±, the ground state 
consists of spin singlet dimers on each rung, as in the 
Heisenberg ladder. The lowest spin excitation is a triplet 
excitation on one rung, propagating via the coupling be- 
tween the rungs. There are also hole excitations with 
spin-i, which carry both spin and charge. The lowest 
2-hole excitation consists of a singlet hole pair on one 
rung, which again develops into a band by propagation 
along the ladder. Relative to the 2-hole ground state, this 
then corresponds to a gapless band of particle- hole charge 
excitations: thus the system is said to be in a Luther- 
Emery C1S0 phase (1 gapless charge mode, gapless 
spin modes). The spin gap evolves discontinuously away 
from half-fillingjll, because there is a triplet particle-hole 
excitation of lower energy than the triplet magnon state 
referred to above. 

A phase diagram in the J/t versus electron density n 
plane, for the case of isotropic couplings, has been pro- 
posed by Poilblanc et alH and by Miiller and RiccM 
For J/t > 2, phase separation is predicted to occur0. 
For J/t < 2 and low to moderate doping the system is in 
a C1S0 phase, and crosses to a C1S1 phase (i.e. the spin 
gap vanishes) for higher doping. For small J/t and low 
doping Miiller and Ricc0 also find ferromagnetic (Na- 
gaoka) and C2S2 phases. Superconducting pairing corre- 
lations occur away from half-firiingQ'BH 

Jurecka and BrenigS have recently constructed an an- 
alytic model of the 2-hole bound states on the t — J lad- 
der, starting from a rung dimer ground state as discussed 
above, and including both 1-hole pseudofermion elemen- 
tary excitations and 2-holc singlet excitations on a single 
rung. They find a low-lying S — bound state in good 
agreement with exact diagonalization calculations, but 



the triplet gap is underestimated as the isotropic limit is 
approached. 

The outline of the paper as follows: in Sect. II we 
briefly describe our series expansion method, while Sec- 
tion III present results for pairs of 'even parity' and 'odd 
parity' excitations respectively, A rich and complex spec- 
trum of bound states is revealed. Sect. IV discusses 4- 
hole bound states formed from two hole-pairs, each on a 
single rung, and this allows some discussion of phase sep- 
aration. Sect. V give a summary and discussion. All the 
bound states discussed in this paper involve identical 'sin- 
gle particle' components; bound states of non-identical 
components will be treated in a separate paper. 

II. METHOD 

We employ a "rung basis" , in which the 2nd and 4th 
terms in (|l|) form the unperturbed Hamiltonian Hq, and 
the remaining terms are treated perturbatively. That is, 
the Hamiltonian (|l|) is rewritten as 



H = H + xV 



(2) 



where 



-y^ P ^li,a^2,a + H.c.)P} (3) 

i . a 



V — Jj_[y^(S^q • Sj + i iQ — -n^ a n^ + i ;Q .) 

i,a 

~V2 E P (4,a,**+l,a,* + H.C.)P] (4) 

z,a,(X 

and x = J/J±, yi = t±/J±, and 7/2 = t/J. 

The eigenstates of Hq are direct products constructed 
from the nine possible rung states, which are listed in Ta- 
ble I for easy reference. In the half-filled case, the lowest 
energy rung state is a spin-singlet state (|tl)— UT) VV2- 
The ground state of the unperturbed Hamiltonian, at 
half-filling, is then a direct product of these spin-singlet 



3 



states on each rung. The lowest spin excitation, at half- 
filling, consists of a spin-triplet excitation on one rung, 
which propagates coherently along the ladder. One and 
two-hole charge excitations are created by removing one 
or two electrons from a rung state and allowing these 
to propagate, via the i-hopping term. In our previous 
workEl, we have studied the 1-hole properties, and the 
properties of 2-holes sitting on the same rung. The study 
of the dynamic properties of 2-holes sitting on different 
rungs is much more complicated: these have been stud- 
ied by Jurecka and BrenigS and by Troyer, Tsunetsugu 
and RiccEl. 

Recently, we have developed strong-coupling series 
expansion methods to study two-particle spectra of 
quantum lattice modelsEl. These methods allow us 
to precisely determine the low-lying excitation spec- 
tra of the models at hand, including all two-particle 
bound/antibound states, at least for small x. In this pa- 
per, we apply this new technique to study the dynamic 
properties of 2-holes in the t — J ladder. 

As discussed in our previous paperEl, to compute the 
two-particle properties, we first calculate an effective 
Hamiltonian in the two-particle sector 

S 2 (i,j;k,l) = (k,l| J ff off |i,j) , (5) 

and then calculate the irreducible two-particle matrix el- 
ement 

A 2 (i,j;k,l) = ^(i,j;k,l)-^ B (*i, k *j,i + *i,i<S llk ) 
-Ai(i,k)*j,i- Axftl)^ 
-Aidk^u-AiCU)^, (6) 

where 5 refers to a Kronecker delta function and Ai is 
the one-particle irreducible matrix clement. Once the ef- 
fective two-particle Hamiltonian is known, we still have 
to solve the Schrodinger equation. The two particle con- 
tinuum is delimited by the maximum (minimum) energy 
of two single particle excitations whose combined mo- 
mentum is the center of mass momentum. There may be 
multiple solutions to the Schrodinger equation above or 



below the two-particle continuum. Those solutions with 
energy below the bottom edge of the continuum are the 
bound states, while those with energy higher than the 
upper edge of the continuum are the antibound states. 
The binding energy is defined as the energy difference 
between the lower edge of the continuum and the energy 
of the bound state, while the antibinding energy is de- 
fined as the energy difference between the upper edge of 
continuum and the energy of an antibound state. 

To compute the perturbation series we choose values 
for the free parameters yi , ?/2 , and derive expansions in 
powers of x. Series have been computed to order x . 
The calculations involve (trivial) 1-dimensional clusters 
up to 12 rungs, and are limited to this order by computer 
memory constraints. To avoid excessive data here we 
mainly consider the case that y\ = yi- The series for 
the energies of the lowest singlet bound states S\ and Si 
and triplet bound state T\ (and also the lower edge of 
the continuum) at band maximum k = ir for yi = y% = 
0.1 are given in Table ||. Other series are available on 
request. 

In the small x limit our strong coupling series are 
highly accurate and we can find all details of various 
two-particle bound (and antibound) states. The over- 
all 2-particle spectrum is much richer than that ob- 
tained in previous studies. Several singlet and triplet 
bound/antibound states are found. The number of 
bound states depends on the coupling constants as well 
as the wavevector. 

Previous numerical studieJll'S on small lattices have 
concentrated on the coupling y\ = 1/3, yi = 10/3, and 
t = tx, and found one singlet bound state in the even 
parity channel. Unfortunately, our series expansions do 
not converge well at this coupling, so we can not make 
a direct comparison. In order to check the correctness 
of our results, we also did some finite lattice calculations 
on small systems, and find very good agreement for the 
lowest energy bound state. The finite lattice calculations 
on small systems certainly are not accurate enough, how- 
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ever, to see other bound states sitting very close to the 
lower edge of the continuum. 

III. TWO-HOLE BOUND STATES 



A. Even-parity channel 




0.2 0.4 0.6 



J/J, 

FIG. 1: A diagram for the number of singlet bound states 
for the even-parity channel in the plane of t/J(— t±/J±) and 
J/J±. 




FIG. 2: A diagram for the number of triplet bound states 
for the even-parity channel in the plane of t/J(— t±/J±) and 
J/J±. 

Among the nine eigenstates for the unperturbed 
Hamiltonian H on each rung, there are the following 
two spin-i electron- hole bonding states^ 

4(|0i)+|4 0» 



(7) 

^(|0T)+|T0», 

which have 'even parity' with respect to reflection across 
the ladder. The eigenenergy for these two states is — 1±. 
In this section, we discuss the bound states formed from 
two of these excitations. 

We find that the number of bound states and an- 
tibound states depends strongly on J/J±_ and t/J(= 
t±/J±). Ignoring the antibound states, Fig. |TJ show the 
number of singlet (S — 0) bound states in the plane of 
t/J(= t±/J±) and J/Jj_. We can see that for large t/J 
(t/J ~ 0-5) there is only one singlet bound state, while 
around t/J ~ 0.3, there is a region where there are two 
singlet bound states, and for smaller t/J and J / J±_ > 0, 
there is a region where there are more than 2 singlet 
bound states. In this region, the number of bound state 
increases as J/J± increases. Note that all these are in 
addition to the singlet state consisting of a pair of holes 
on the same rung. 

Fig. H shows the number of triplet (S = 1) bound 
states in the same plane. We can see that for large t/J 
(t/J > 0.5) there are no triplet bound states, whereas 
around t/ J ~ 0.2 and small J/J±, there is a region where 
there is only one triplet bound state, and for smaller t/J 
and J/ J±_ > 0, there is a region where there are multi- 
ple triplet bound states. The existence of these multiple 
bound states has not been noted before, to our knowl- 
edge. Note also that the triplet bound states appear to 
vanish as the isotropic case J/Jj_ is approached. 

Now let us look at a particular case t/J = t±/J± = 0.1 
(i.e. small hopping parameters) and J/Jj_ = 0.2. The 
two-particle excitation spectrum is shown in Fig. |^. One 
can see from these graphs that there are three singlet 
bound states (Si, S2 and S3) and one triplet bound state 
(Ti) below the continuum, and two triplet antibound 
states (T2 and T3) above the continuum, in addition to 
the state So corresponding to a hole pair on a single 
rung. All the dispersion bands are very flat, as expected 
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FIG. 3: The excitation spectrum for the even parity chan- 
nel with t/J = t±/J± =0.1 and J/J± = 0.2. Beside the 
two-particle continuum (gray shaded), there are three singlet 
bound states (Si, S2 and S3) and one triplet bound state (Ti) 
below the continuum, and two triplet antibound states (T2 
and T3) above the continuum. The insets enlarge the region 
near k = ir so we can see S2, S3, Ti and T2, but the antibound 
state T3 is still indistinguishable from the continuum. Curve 
A is the dispersion of the one-hole bonding state, while curve 
So is the dispersion of two-holes sitting on the same rung. 
The cross points are the results of exact diagonalization for 
the 8-rung finite lattice with periodic boundary condition. 

for small hopping parameters. The singlet bound state 
S\ exists for the whole range of momenta (its coherence 
length L is finite also for the whole range of momenta), 
while other bound/antibound states exist only in a lim- 
ited range of momenta near k = ir. Figure |^ shows the 
inverse of the coherence length 1/L versus momentum k 
for these various bound states. At k = w, the coherence 
length L for S\ and T\ is 1, for S2 and T2, it is about 2, 
S3 has coherence length L about 3, while for T3 the co- 
herence length L is very large, about 20 (its antibinding 
energy at k = ir is very small, E\, /J± = 8.9 x 10 -8 , so T3 is 
indistinguishable from the two-particle continuum in Fig. 
|). This means that the formation of bound/antibound 
states Si, Ti , S%, T2, S3 is largely due to the interac- 
tion of two bonding rungs separated by distances 1, 1, 2, 
2, and 3 respectively, while the formation of T3 is much 
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FIG. 4: The inverse of the coherence length 1/L versus mo- 
mentum k for three singlets (Si, S2 and S3), one triplet (Ti) 
bound state, and two triplet (T2 and T3) antibound states 
for the even parity channel at t/J — t±/J± — 0.1 and 
J/ J± — 0.2. The inset enlarges the region near k = n. 

more complicated. For comparison, also shown in Fig. |3| 
is the dispersion for the one-hole bonding state, and the 
dispersion of two-holes sitting on the same rungH: here 
two holes sitting on the same rung have lower energy. 

To check the correctness of our series calculations, we 
also performed an exact diagonalization for a finite lattice 
of 8 rungs with periodic boundary condition for this set 
of parameters. The results for k = 0, ir are also shown in 
Fig. ||. There is remarkable agreement for the state Si, 
but the finite lattice calculation is not accurate enough 
to give the other bound states sitting very close to the 
continuum. 

We now study the behaviour as a function of the cou- 
pling ratio J/Jxj fixing t/J = tx/J± = 0.1 and k = ir. 
The binding/antibinding energies Ef, and the inverse of 
the coherence length 1/L versus J/J± are shown in Figs. 
|B| and|^. In the small J limit, there are two singlet bound 
states Si and £2 (with coherence lengths 1 and 2, re- 
spectively), one triplet bound state Ti (coherence length 

1) , and one triplet antibound state Ti (coherence length 

2) . As J increases, more and more singlet and triplet 
bound states appear: we find that the singlet /triplet 
bound/antibound states S3, S4, S§ and T3 appear at 
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J/Jj_ about 0.0480, 0.39, 0.55 and 0.180 respectively. To 
find out the behaviour near the threshold, we plot in Fig. 
the square root of the binding energy (Et,/J±) 1 / 2 and 
the inverse of the coherence length 1/L versus J/J± for 
singlet bound state S3. We can see that near threshold, 
the curves are nearly straight lines, which implies that 
near threshold Eb oc (x — x c ) 2 and L oc (x — x c )~ l , cor- 
responding to "threshold indices" 2 and -1, respectively, 
where x = J/J±. For a fixed value of J/J±, the threshold 
behaviour as a function of momentum is Eb oc (k — fc c ) 2 , 
L oc (k — k c )~ l . We expect similar behaviour for other 
bound/antibound states (but see Sect. IV). 
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FIG. 5: The binding energy Eb and the inverse of the coher- 
ence length 1/L versus J/J± for the singlet bound states in 
the even parity channel at t/J = t±/J± — 0.1 and k = tt. 
For the binding energy of Si and 5*2, the results of several 
different integrated differential approximants to the series are 
shown, while for other curves, the results of the three highest 
orders are plotted. 

Next we consider the case t/ J = t±/ J± = 1 (stronger 
hopping terms) and J/J± = 0.2. The two-particle exci- 
tation spectrum and the inverse of the coherence length 
1/L versus momentum k are shown in Figs. ^ and |^. 
Here there is only one singlet bound state Si below the 
continuum, and one triplet antibound state (Ti) above it. 
Both 5*1 and T± exist only in a limited range of momenta 
near k = it, and their coherence lengths L at k — n are 
about 1. As comparison, also shown in Fig. @ is the dis- 
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FIG. 6: The binding/antibinding energy Eb and the in- 
verse of the coherence length 1/L versus J/J± for triplet 
bound/antibound states in the even parity channel at t/J = 
t±/J± = 0.1 and k — tt. For the binding/antibinding energy 
of Ti and T2, the results of several different integrated differ- 
ential approximants to the series are shown, while for other 
curves, the results of the three highest orders are plotted. 
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FIG. 7: The square root of the binding energy (Eb/ J±) 1 ^ 2 
and the inverse of the coherence length 1/T versus J/J± for 
the singlet bound state S3 in the even parity channel at t/J = 
t±/J± = 0.1 and k — tt. 

persion for the one-hole bonding state, and the dispersion 
of two- holes sitting on the same runglEI, So- In this case 
the state So actually lies above the state Si, and even 
above the 2-particle continuum. 

If we fix the values t/J = t±/J± = 1 and fc = 7r, 
the binding/antibinding energies Ej, and the inverse of 
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FIG. 8: The excitation spectrum for the even parity channel 
with t/J — t±/J± — 1 and J/J± = 0.2. Beside the two- 
particle continuum (gray shaded), there is one singlet bound 
state (Si) below the continuum, and one triplet antibound 
state (Ti) above the continuum. The insets enlarges the re- 
gion near k = n to show Ti above the continuum. The curve 
labeled A is the dispersion of one hole, while the curve labeled 
So is the dispersion of two holes on the same rung. 




k/7T 

FIG. 9: The inverse of the coherence length 1/L versus mo- 
mentum k for the singlet Si and one triplet antibound state 
7i in the even parity channel at t/J — t±/J± = 1 and 
J/Jx = 0.2. 

the coherence length 1 /L versus J/ Jj_ are shown in Fig. 
|Io| . One can see that the singlet bound state Si and 
the triplet antibound state Ti exist in the small J limit, 
but as J increases, no more new bound/antibound states 
appear. 
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T, (n=2, antibound) 

H 1 1 1 1 1 



even parity k-TT 



t/J=t i /J 1 =l 
even parity k = 7T 



J/J 1 

FIG. 10: The (anti)-binding energy Eh and the inverse of the 
coherence length 1/L versus J/J± for singlet /triplet (anti- 
bound states in the even parity channel with t/J = t± /J± — 1 
and k = n. For the (anti)-binding energy, the results of several 
different integrated differential approximants to the series are 
shown, while for the inverse of the coherence length 1/L, the 
results of the three highest orders are plotted. 

Sometimes a bound state can appear in the regions 
near both k = and k — n, but not at intermediate 
momenta. For example for t/J = t±/J± = 0.425, the 
binding energy and coherence length L for J/J± = 0.025, 
0.05 and 0.075 are shown in Fig. [□]. For very small J j J± 
(0.025 and 0.05 in the figure), the singlet bound state 
exists over the whole range of momenta, but for larger 
J/J± (0.075 in the figure), this singlet bound state only 
exists in the regions near both k = and k = n, but not 
in the intermediate regions of momentum. 

As already seen above, the lowest energy 2-hole bound 
state may correspond to a pair of holes on the same rung 
(Sq), or two holes on different rungs (Si), depending on 
the parameters t/J(= t±/J±) and J/J±. The boundary 
between these two situations is shown in Fig. [l^: we will 
refer back to this diagram later on. 

B. Odd-parity channel 

Among nine rung eigenstates for the unperturbed 
Hamiltonian Hq in each rung, there are the following two 
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J/J =0.026 \. 




FIG. 11: The binding energy E\> and the inverse of the co- 
herence length 1/L versus k for the first singlet bound state 
in the even parity channel at t/J — t±/J± = 0.425 and 
J/J± = 0.025 (solid lines), 0.05 (dotted lines), and 0.07 
(dashed lines). The results of the three highest orders are 
plotted. Note that for these couplings there is a second sin- 
glet bound state and a triplet bound state near k = tt, which 
are not plotted. 
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FIG. 12: Boundary between different configurations of the 
lowest-lying 2-hole bound state. 

spin-i odd-parity electron-hole antibonding states 
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(8) 



The eigen energy for these two states is t±. In this sec- 
tion, we discuss bound states of these two 1-particle ex- 
citations. 



one S-0 bound state 




two S = bound states 



>2 S = bound states 



0.1 0.2 0.3 

J/J, 

FIG. 13: The number of singlet bound states in the odd-parity 
channel in the plane of t/J(= t±/J±) and J/J±. 



odd parity 



one S-l bound state 



no S=l bound states 




FIG. 14: The number of triplet bound states in the odd-parity 
channel in the plane of t/J(= t±/J±) and J/J±. 

As for the case of even parity, the number of bound 
states and antibound states depends on J/J± and tj J(= 
tj_/Jj_). Fig. |l3| shows the diagram for the number of 
singlet (S = 0) bound states, where we can see that for 
large t/J (t/J > 0.5) there is only one singlet bound 
state, while around t/J ~ 0.3, there is a region where 
there are two singlet bound states. For smaller t/J and 
J / 'J± > 0, there is a region where there are more than 
two singlet bound states, and in this region, the number 
of bound states increases as J / J± increases. 

Fig. [li] shows the diagram for the number of triplet 
(S = 1) bound states. For t/J> 0.5, there is one triplet 
bound state, while elsewhere no triplet bound state ex- 
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ists, except in a region at very small t/J. 
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FIG. 15: The excitation spectrum in the odd parity channel at 
t/J = t±/J± = 0.1 and J/Jx = 0.2. Beside the two-particle 
continuum (gray shaded) , there are three singlet bound states 
(Si, S2 and S3) below the continuum, and three triplet an- 
tibound states (Ti, T2 and T3) above the continuum. Two 
lower insets enlarge the region near k = n to show S2, S3, Ti 
and T2 below/above the continuum, but the antibound state 
T3 is still indistinguishable. To see T3, we plot the antibinding 
energy versus k/n, as the top inset. 

If we take the values t/J = t±/J± = 0.1 and J/J± = 
0.2, the two-particle excitation spectrum and the in- 
verse of the coherence length 1/L versus momentum k 



are shown in Figs. |15j and 16. There are three singlet 
bound states (Si, S2 and 53) below the continuum, and 
three triplet antibound states (Ti, T2 and T3) above the 
continuum. The singlet bound state Si exists for the 
whole range of momenta (its coherence length L is fi- 
nite also for the whole range of momenta), while other 
bound/antibound states exist only in a limited range of 
momenta near k — n. At k = tt, the coherence length 
L for Si and Ti is 1, for S2 and T 2 it is about 2, while 
S3 and T3 have coherence length L about 3. This means 
that the formation of bound/antibound states Si, Ti , 
S2, T2, S3 and X3 is largely due to the interaction of two 
antibonding rungs separated by distances 1, 1, 2, 2, 3, 
and 3 respectively. 

If we fix the values t/J = tx/Jx = 0.1 and k = tt, 




FIG. 16: The inverse of the coherence length 1/L versus mo- 
mentum k for three singlet (Si, S2 and S3), and three triplet 
(Ti, T2 and T3) antibound states for the odd parity channel 
at t/J — tx/Jx = 0.1 and J/Jx = 0.2. The inset enlarges 
the region near k — tt. 

the binding/ antibinding energies Eb and the inverse of 
the coherence length 1/L versus J/J± are shown in Figs. 
and H One can see that at the small J limit, we 
have two singlet bound states Si and S2 (with coherence 
lengths 1 and 2, respectively), and two triplet antibound 
states Ti and T2 (with coherence lengths 1 and 2, respec- 
tively). As J increases, more and more singlet bound 
state and triplet antibound states turn up, and in our 
calculations, we find that states S3, S4, and T3 appear at 
J/J± about 0.05, 0.12, and 0.07 respectively, while S2, 
S3, S4, and Ti disappear at J/J± about 0.32, 0.29, 0.13 
and 0.46, respectively. Thus state S4 pops out below the 
continuum only very briefly. The threshold index for the 
binding energy of S2 and Ti is about 1, rather than 2, as 
found by Dlog Pade approximants to the series. 

For larger hopping terms t/J = tx/Jx = 1 and 
J I Jx = 0.2, the two-particle excitation spectrum and 
the inverse of the coherence length 1/L versus momen- 
tum k are shown in Fig. [l9|. Now there is one singlet 
bound state Si and one triplet bound state (Ti) below 
the continuum, but no antibound state. Both Si and Ti 
exist only in a limited range of momenta near k = w, and 
their coherence lengths L at k = tt are about 1. 




o.a - t/j-tyj^o.i 

: odd parity S=0 k=7T 




0.1 0.2 0.3 



J/J ± 

FIG. 17: The binding energy Et and the inverse of the coher- 
ence length 1/L versus J/J± for singlet bound states in the 
odd parity channel at t/J = t±/J± — 0.1 and k = n. For the 
binding energy of Si and 5*2, the results of several different 
integrated differential approximants to the series are shown, 
while for other curves, the results of the three highest orders 
are plotted. 

If we fix the values t/J = tj_/J± = 1 and k = tt, 
the binding energies and the inverse of the coherence 
length 1/L versus J / J± are shown in Fig. [2(]. The singlet 
bound state Si and the triplet bound state Ti exist as 
long as J =/= , but no more new bound/antibound states 
appear as J increases. 

IV. FOUR-HOLE BOUND STATES AND PHASE 
SEPARATION 

Among nine rung eigenstates for the unperturbed 
Hamiltonian Hq on each rung, there is one state of a two- 
hole pair with eigenenergy 0. In this section, we discuss 
bound states of 2 two- hole pair excitations, correspond- 
ing to the half filled system doped with 4 holes. 

As for the case of 2-hole bound states, the number 
of bound states and antibound states depends on J / J±_ 
and t/J(= t±/J±). Fig. pl| shows the number of singlet 
(S = 0) bound states in various regions of parameter 
space. For large t/J (t/J > 0.3) and large J/J±, there 
is only one singlet bound state, while elsewhere there 
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0.1 0.2 0.3 0.4 

J/J ± 



FIG. 18: The antibinding energy Et, and the inverse of the 
coherence length 1 / L versus J / J± for triplet antibound states 
in the odd parity channel with t/J = t±/J± = 0.1 and k — tt. 
For the antibinding energy of Ti and T2, the results of several 
different integrated differential approximants to the series are 
shown, while for other curves, the results of the three highest 
orders are plotted. 

are two singlet bound states, except for a peculiar very 
narrow region (given by the bold solid line) where there 
is only one bound state, and a small region at low t/J 
(labelled by > 2) where there are more than two bound 
states. 

If we take the values t/J = t±/J± = 0.1 and J/J± = 
0.15 (we take J/J± = 0.15, rather than J/J± = 0.2 as 
before, because at J/J± = 0.2 we have only one bound 
state), the two-particle excitation spectrum and the in- 
verse of the coherence length 1/L versus momentum k 
are shown in Figs. ^ and |2^. There are two singlet 
bound states (Si and S2) below the continuum. The 
singlet bound state Si exists for the whole range of mo- 
menta (its coherence length L is finite also for the whole 
range of momenta), while bound state S2 exists only in 
a limited range of momenta near k = tt. The coherence 
length L for Si is 1, while the coherence length L for S2 
at k = n is about 2. This means that the formation of 
bound states Si and S2 is largely due to the interaction 
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I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 
t/j=t J /j 1 =i. j/j^o.a 

..... odd parity 




FIG. 19: The excitation spectrum and the inverse of the 
coherence length 1/L in the odd parity channel at t/J = 
t±/J± = 1 and J/J± = 0.2. Beside the two-particle contin- 
uum (gray shaded), there are one singlet (Si) and one triplet 
bound state (Ti ) below the continuum. 




J/J, 

FIG. 20: The binding energy Et and the inverse of the co- 
herence length 1 / L versus J / J± for singlet and triplet bound 
state in the odd parity channel at t/J — t±/J± = 1 and 
k = 7r. For the binding energy, the results of several different 
integrated differential approximants to the series are shown, 
while for the coherence length, the results of the three highest 
orders are plotted. 

of two hole-pair rungs separated by and 1 singlet rung, 
respectively. 

If we fix the values t/J = t±/J± = 0.1 and k — ir, the 
binding/antibinding energies Ej, and the inverse of the 
coherence length 1/L versus J/J± are shown in Fig. |24|. 
One can see that in the small J limit, we have two sin- 



0.4 



two bound states 



one bound state 



two bound states 




FIG. 21: The number of singlet 4-hole bound states in the 
plane of t/J(— t±/J±) and J/J±. Note that the bold solid 
line is the region where there is only one bound state, whereas 
> 2 is the region that there are more than 2 bound states. 

2.34 



2.32 



< 



0.2 0.4 0.6 0.8 1 

k/7T 

FIG. 22: The excitation spectrum for 4-hole bound states at 
t/J = t±/J± = 0.1 and J/J± = 0.15. Beside the two-particle 
continuum (gray shaded), there are two singlet bound states 
(Si and S2) below the continuum. The inset enlarges the 
region near k — tv to show S2 below the continuum. 

glet bound states S\ and S2 (with coherence length being 
1 and 2, respectively). As J increases, S2 disappears at 
J I 'J± ~ 0.193, while a new antibound state S3 and a new 
bound state £4 appear at J/J± ~ 0.205 and 0.235, (co- 
herence length at larger J being 2 and 3), respectively. 
This interval corresponds to the peculiar region given by 
the bold line in Fig. |l|. The threshold index as the 
binding energy of disappears appears to be about 1 , 
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FIG. 23: The inverse of the coherence length 1/L versus mo- 
mentum k for two singlet (Si and S2) 4-hole bound states at 
t/J — t±/J± = 0.1 and J/J± — 0.15. The inset enlarges the 
region near k — it. 

rather than 2. Estimates for the threshold point and in- 
dex from the [n/m] Dlog Padc approximants to the series 
for the binding energy are given in Table 



III 



We cannot 



calculate a series directly for the coherence length, but 
from our numerical calculations, the threshold behaviour 
for the coherence length of S% is also unlikely to be the 
standard L cx (x — x c )~ l , but seems to be a weaker diver- 
gence. The dynamical reason for this unusual threshold 
behaviour is unclearEl. Note that even in this case, the 
threshold behaviour as a function of momentum still re- 
tains the form: Ef, cx (k — k c ) 2 , L cx (k — fc c ) _1 . 

For larger hopping terms t/J = tj_/Jj_ = 1 and 
J/J±_ — 0.15, the two-particle excitation spectrum and 
the inverse of the coherence length 1/L versus momentum 



k are shown in Figs. 25 and E6L Again there are two sin- 



glet bound states £1 and S2 below the continuum, but 
no antibound state. The singlet bound state Si exists 
for the whole range of momenta (its coherence length L 
is about 1 also for the whole range of momenta), while 
bound state S2 exists only in a limited range of momenta 
near k = it (its coherence length L is about 2 at k = tt). 

If we fix the values t/J = t±/J± = 1 and k = 7r, 
the binding energies Ey, and the inverse of the coherence 
length 1/L versus J/J± are shown in Fig. E7| The two 
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1 f S ° 



0.0001 

8x10"* 
^ 6xl0~ 5 
H-4X10- 

2xl0" 6 


1 

0.8 
0.6 
0.4 
0.2 


0.1 0.2 0.3 0.4 0.5 

J/J, 

FIG. 24: The binding energy Eb and the inverse of 
the coherence length 1/L versus J/J± for singlet 4-hole 
bound/ antibound states at t/J = t±/J± = 0.1 and k — n. 
For the binding energy of Si and 52, the results of several 
different integrated differential approximants to the series are 
shown, while for other curves, the results of the three highest 
orders are plotted. 



< 



2.382 
2.38 
2.378 



2.376 



0.9 0.92 0.94 0.96 0.98 1 



0.2 0.4 0.6 0.8 1 

k/7T 

FIG. 25: The excitation spectrum for 4-holes bound state at 
t/J — t±/J± — 1 and J/J± = 0.15. Beside the two-particle 
continuum (gray shaded), there are two singlet (Si and S2) 
below the continuum. Curve A is the dispersion of the two- 
holes sitting in the same rung. 



singlet bound states 5*1 and S% exist as long as J / , 
but as J increases, no more new bound/antibound states 
appear, but S2 disappears at J / 'J± ~ 0.22. 
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(L/2)[A£ (4 holes) - 2A£ (2 holes)] 



(10) 



0.2 - 



0.2 0.4 0.6 O.f 



k/TT 

FIG. 26: The inverse of the coherence length 1/L versus mo- 
mentum k for two singlet (Si and S2) 4-holes bound states at 
t/J = t±/J± = 1 and J/J± = 0.15. The inset enlarges the 
region near k — ty. 




FIG. 27: The binding energy E\> and the inverse of the co- 
herence length 1 / L versus J / J± for two singlet 4-holes bound 
states at tj J — t±/ J± = 1 and k = tt. For the binding energy, 
the results of several different integrated differential approxi- 
mants to the series are shown, while for the coherence length, 
the results of the three highest orders are plotted. 



The compressibility k is given by 



,d 2 



dp 2 



(9) 



where e(p) is the energy density of the ladder with hole 
density p. For a finite system at half-filling this can be 
replaced by the discrete version 

k.- 1 = (L/2)[E (- filled) - 2E Q (2 holes) + E (4 holes)] 



where for a ladder of L rungs, AEq is the minimum en- 
ergy gap for 2 or 4 holes, and this minimum energy gap 
is at momentum k = 0. Now if the 2-hole pair on a sin- 
gle rung is the lowest-energy state in the 2-hole sector, 
and at momentum k = two 2-hole rungs form a bound 
state in the 4-hole sector, then it immediately follows 
that kT 1 as defined by Eq. [l^ diverges to negative infin- 
ity as L — * 00, signalling phase separation. Comparing 
Figures [2l] and [l^, we see that these conditions apply, 
and phase separation occurs, over the whole of the lower 
region of Figure [l2[ 

Outside this region, we know the energy of the lowest 
state in the 2-hole sector, consisting of a bound pair of 
holes on different rungs (denoted S%, previously). The 
bound state of two 2-hole rungs may not be the lowest- 
energy state in the 4-hole sector; but at least it allows us 
to place an upper limit on kT 1 as defined by Eq. [l^. We 
find that this limit is negative, and therefore phase sep- 
aration must occur, for t/J(= Jj_) anywhere below a 
value about 0.5, independent of J / J±_ - although due to 
numerical uncertainty, this estimate cannot be made very 
precise. This is in good agreement with previous esti- 
mates: Rommer et a/@, for instance, find the phase sep- 
aration boundary at half-filling to lie at J/t = 2.156(2), 
or t/J — 0.464, from a density matrix renormalization 
group calculation for the isotropic ladder. 

V. DISCUSSION AND SUMMARY 

We have used linked-cluster series expansion methodsi 
to study 2-hole bound states in the t — J ladder. The se- 
ries only converge well for J/J± ^ 0.5, so we are unable 
to say very much about the isotropic case J = J±. Dis- 
persion relations and coherence lengths have been com- 
puted for various bound states in even- and odd- parity 
channels. A rich spectrum of these states has emerged. 

The most striking feature of the results is the appear- 
ance of multiple bound states in both singlet S = 
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and triplet S = 1 channels as the ratio t/J(= t±/J±) 
goes to zero for finite J/J±. What is the cause of this 
phenomenon? The limit t/J — > corresponds to the 
antifcrromagnetic Heisenberg ladder with 2 static holes. 
Beyond states Sq and Si, all the multiple bound states 
(S2, ■ • •) only emerge near k = n, and all have small bind- 
ing energies: their dynamical significance is unclear. 

Another interesting feature is that in the even-parity 
channel, at least, all the triplet 5 = 1 bound states 
seem to disappear as the isotropic limit J / J±_ — > 1 is 
approached, leaving only singlet 5 = bound states. 
Other studiesil have found the same thing. We have no 
explanation for this phenomenon, either. 

Jurecka and BrenigS have put forward a semi-analytic 
theory of 2-hole excitations in the t — J ladder, based 
on the dimerized rung picture, with elementary excita- 
tions consisting of the single-hole pseudofcrmion excita- 
tions and the 2-hole singlet excitation on a single rung, 
plus singlet and triplet "bond boson" operators for spin 
excitations, and applying a linearized Holstein-Primakoff 
approach. They present a spectrum showing only one 
5 = and one 5=1 bound state, at couplings which 
are not accessible to our series approach. 

The question of phase separation has also been ex- 
plored. We have been able to calculate the energy of 
a 4-hole state consisting of two 2-hole rungs, giving an 
upper limit on the inverse compressibility at half-filling. 
Hence it can be shown that phase separation occurs for 
t/J < 0.5, in agreement with other studies0. Most of 
the interesting spectral features occur, unfortunately, in 
the phase separated region. 

We have also explored the threshold behaviour of these 
bound/antibound states as they emerge from the contin- 
uum. As a function of momentum, the binding energy 
near threshold behaves as Ef, cx (k — k c ) 2 , and the co- 
herence length as L cx (k — fc c ) _1 in every case we have 
studied so far. This appears to be standard: for instance, 
a model of the bound states of two magnons by Sushkov 
and Kotov0 gives this behaviour. As a function of cou- 



pling x = J/J±, the typical threshold behaviour is sim- 
ilar, Eb cx (x — x c ) 2 , L cx (x — x c ) _1 . We have found 
examples, however, where the binding energy appears to 
behave as E^ cx (x — x c ), see Sect. IIIB and IV. 

In conclusion, then, we have discovered a rich and di- 
verse spectrum of 2-hole bound states in the t — J ladder 
at half-filling. They exhibit a number of interesting fea- 
ture which call for a theoretical explanation; and fitting 
these data will provide a stringent test of any theoretical 
model of the system. It would also be interesting if it 
were possible to test these features experimentally. 
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TABLE I: The nine rung states and their energies at t = J = 0. 



INo. 


Eigenstate 


Eigenvalue 


Name 


1 


75 (IU>- lit)) 


T 


singlet 


2 


ii) 





triplet (Sjot = -1) 


3 


75 (lti)+ lit)) 





triplet (S tot = 0) 


4 


Itt) 





triplet (6 tot = 1) 


o 


l nn\ 
uu; 


n 

u 


iiuie-pdii bingiet 


6 


75(|0i)+|i0» 


-t± 


electron-hole bonding (S^ot = — 5) 


7 


75(|0t)+|t0)) 


-t± 


electron-hole bonding (£tot = 5) 


8 


7^(1 1)- Ii0» 


t± 


electron-hole antibonding (Stot = — |) 


9 


75(1 °t)- IT o» 


t± 


electron-hole antibonding (5 t z ot = |) 



TABLE II: Series coefficients for dimer expansions for the energy gap E/J± of two singlet bound states (Si and 52), one triplet 
bound state (Ti), and the lower edge of the continuum (C;) at k = it, for the bound states of two even parity holes of the t — J 
ladder with y 1 = y 2 = 0.1. Nonzero coefficients x n up to order n = 11 are listed. 



n 


E S JJ± for Si 


E S JJ± for S 2 


E T JJ± for Ti 


E Cl /J± for Ci 





1.8000000000 


1.8000000000 


1.8000000000 


1.8000000000 


1 


5.0000000000X10" 1 


1.0000000000 


1.0000000000 


1.0000000000 


2 


2.2500000000xl0 _1 


8.5000000000X10" 1 


8.0000000000X10" 1 


8.5000000000X10" 1 


3 


4.2708333333X10" 1 


3.2811666667X10" 1 


3.8020833333xl0 _1 


4.5625000000X10" 1 


4 


9.8972656250X10" 1 


-5.0189024889xl0 _1 


-2.1309606481xl0 _1 


-2.0194444444xl0 _1 


5 


-1.4465469473 xlO" 1 


-1.0020098268 


-6.5367901837xl0 _1 


-7.4625737847xl0~ 1 


6 


-3.6759280211 


-5.2752447741 xlO -1 


-5.1516620601xl0 _1 


-5.7308628683X10" 1 


7 


-1.8432595464 


9.3228166551 xlO -1 


2.9350608457xl0 _1 


4.7837543431 xlO -1 


8 


1.2540109894X10 1 


2.4273057426 


1.1202826230 


1.5454790111 


9 


1.1255237404X10 1 


2.2797858202 


8.6634384749xl0 _1 


1.0984536226 


10 


-4.7671925391 xlO 1 


-5.2477491651 xlO -1 


-8.6710572784xl0 _1 


-1.4110630505 


11 


-5.6591915055X10 1 


-4.6480530617 


-2.6386765415 


-3.8611918629 



TABLE III: [n/m] D log Pade approximants to the series for bining energy E b /J± at k — ir of singlet bound state S2 of 4-holes 
at t/J = t±/J± = 0.1. An asterisk denotes a defective approximant. 
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[(n - 2)/n] 


[(n - l)/n] 


[n/n] 


[(n + l)/n] 


[(n + 2)/n] 






pole (residue) 


pole (residue) 


pole (residue) 


pole (residue) 


pole (residue) 


n= 


1 






0.19966(1.1704) 


0.19987(1.1741) 


0.19391(1.0401) 


n= 


2 




0.19987 ( 1.1741) 


0.19967(1.1705) 


0.19194(0.9851) 


0.19250(1.0032) 


n= 


3 


0.19415(1.0485)* 


0.19255 ( 1.0055) 


0.19265(1.0088) 


0.19272(1.0117) 


0.19280(1.0154) 


n= 


4 


0.19264(1.0086) 


0.19288 ( 1.0202)* 


0.19296(1.0259)* 


0.19171(1.0808)* 




n= 


5 


0.19295(1.0259)* 


0.19285 ( 1.0185)* 









